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This  paper  studies  the  propagation  of  weakly  nonlinear 
waves.  The  analysis  is  applicable  to  plasma  waves,  surface 
water  waves,  the  interaction  of  laser  beams  with  matter,  particle 
accelerators,  etc.  The  specific  motivation  for  this  work  is, 
however,  the  study  of  internal  waves  in  the  ocean.  Hamilton's 
principle  is  used  to  write  the  fluid  equations  in  Hamiltonian 
form  in  terms  of  linear  eigenmode  amplitudes.  Numerical  studies 
are  made  of  the  effect  of  Fourier  grid  size  and  resonance  widths. 
Statistical  information  is  generated  from  an  ensemble  of  initial 


states  of  the  random  wave  field. 


79  02  09  080 

°7/9?o  _ 


, 


1.  Introduction 

In  this  paper  we  consider  some  techniques  used  to  study  statistical 
properties  of  weak  turbulence  associated  with  nonlinear  wave-wave 
interactions.  Specifically,  we  apply  the  analysis  to  buoyancy  dominated 
turbulence  of  interned  waves  in  the  ocean.  Stable  stratification,  a 
characteristic  of  the  oceans,  implies  an  equilibrium  depth  about  which 
each  fluid  element  oscillates.  The  resulting  "almost  two-dimensional" 
system  avoids  some  of  the  complexity  of  fully  three-dimensional  weak 
turbulence. 

Formally,  our  wave  system  corresponds  to  a set  of  harmonic  oscillators 
with  weakly  nonlinear  couplings.  Similar  dynamical  systems  are  encountered 
in  the  study  of  surface  water  waves  (l),  plasma  waves  (2),  the  interaction 
of  light  with  matter,  particle  accelerators,  etc.  Interest  in  the  pro- 
perties of  oscillators  with  weakly  non-linear  couplings  has  been  stimu- 
lated by  the  KAM  theorem  (3),  the  discovery  of  solitons  (U)  and  other 
studies  (5). 

Relaxation  phenomena  associated  with  linearly  coupled  oscillators 
have  been  studied  by  Mazur  and  Montroll  (6)  and  by  Ford,  Kac  and  Mazur  (7). 
Computations  with  a linear  chain  of  fifty  oscillators  by  Cukier , Shuler , 
and  Weeks  (8)  have  shown  agreement  with  the  predictions  of  a Langevin  model 
for  the  system. 

Weakly  nonlinear  oscillator  systems  have  been  studied  by  several 
statistical  models.  Use  of  random  phase  and  two-time  scale  approximations 
allows  termination  of  a sequence  of  coupled  moment  equations  to  give  a 
Boltzmann-type  transport  equation.  This  approach  has  been  used  by  Hasselman 


j 
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(9)  for  geophysical  phenomena,  Snider  (10)  for  gas  kinetic  theory  and 
Davidson  (2)  for  plasmas.  The  technique  of  Prigogine  and  collaborators 
(ll)  of  simplifying  the  Liouville  equation  by  the  neglect  of  higher  order 
correlations  also  leads  to  the  Boltzmann  transport  equation.  In  a future 
publication  we  will  show  how  the  fluctuation-dissipation  theorem  (12) 
and  the  Krylov-Bogoliubov-Mitropolsky  two-time  perturbation  method  (13,1*0 
can  be  used  to  obtain  Langevin  and  Fokker-Planck  descriptions  of  internal 
wave  turbulence. 

These  models  can  be  compared  with  numerical  solutions  of  the  equations 
of  motion,  providing  evidence  for  the  validity  of  the  approximations  used. 
The  ocean  internal  wave  system  lends  itself  particularly  well  to  numeri- 
cal computation  because  of  the  "almost  two-dimensional"  nature  of  the 
equations.  Recently,  the  Hasselman  transport  theory  has  been  applied 
to  this  system  by  McComas  and  Bretherton  (15)  and  Olbers  (l6).  These 
studies  utilize  the  remarkable  synthesis  of  experimental  observations, 
made  by  Garrett  and  Munk  (17, 18),  who  proposed  an  explicit  equilibrium 
spectrum. 

In  this  paper  an  explicit  Hamiltonian  is  given  which  describes 
the  nonlinear  transfer  of  energy  among  the  linear  eigenmodes  of  the 
internal  wave  field.  A "test  wave"  model  is  developed  which  can  be  used 
to  compare  numerical  results  with  statistical  models.  This  model  describes 
the  propagation  of  a single,  labeled  wave  through  an  ambient  medium. 
Computational  methods  for  integration  of  Hamilton's  equations  are  also 
discussed. 
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2.  Dynamical  Formulation  of  the  Wave  Interactions 


We  consider  as  a model  a "plane”  ocean  of  uniform  depth  and  having 
rectangular  area  I^.  Periodic  boundary  conditions  are  used  at  the  sides, 
and  the  top  and  bottom  surfaces  are  assumed  rigid.  The  quiescent  upper 
surface  coincides  with  the  plane  z = 0 of  a rectangular  coordinate  system. 
The  bottom  is  at  z * -II.  The  equilibrium  fluid  density  is  p(z),  a 
monotonically  increasing  function  of  depth.  The  quantity 


N(z) 


SL  *£. 
p dz 


[2.1] 


is  the  V'aisala,  or  buoyancy  frequency  (g  is  the  acceleration  due  to 
gravity).  The  fluid  is  incompressible  and  inviscid. 

We  introduce  a Lagrangian  and  use  Hamilton's  principle  to  obtain 
the  equations  of  fluid  motion.  The  Lagrange  coordinate  of  a fluid  particle 
at  time  t is  Y(r,t),  where  r * Y(r,0)  is  its  position  at  t = 0.  The 
appropriate  Lagrangian  per  unit  area  is  (20) 

L = / ~ P | Y| 2 + pg-Y  - ~ pf* (YxY)  + P(r)  j^J  J - 1 

[2.2] 

The  first  three  terms  in  the  integrand  represent,  respectively,  kinetic 
energy,  negative  of  gravitational  potential  energy,  and  rotational  energy 
due  to  Corolis  coupling.  The  quantity  "f"  is  twice  the  angular  frequency 
of  the  earth's  rotation.  The  final  term  expresses  the  constraint  due  to 
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incompressibility : P is  a Lagrange  multiplier  and  J is  the  Jacobian 
of  the  transformation  r -+■  Y. 

Hamilton's  principle  states  that  the  functional 


I = / L dt 


[2.3] 


is  stationary  with  respect  to  arbitrary  independent  variations  in  Y and 
P,  which  vanish  at  t^  and  t^  and  at  the  boundaries  of  the  fluid.  In 
particular,  variation  of  P yields  the  incompressibility  condition 


[2.M 


and  variation  of  Y yields  the  equations  of  motion 


pY  - pg  + pf  x Y + — P = 0 


[2.5] 


The  Lagrange  multiplier  P can  therefore  be  identified  with  the  fluid 
pressure.  It  may  be  considered  a function  of  Y and  t in  the  Lagrangian 
since  this  adds  terms  with  zero  variation. 

Rather  than  use  [2.U]  and  [2.5]  directly,  we  develop  a perturbation- 
variation  approach  with  the  aim  of  obtaining  a Hamiltonian  which  describes 
the  lowest  order  nonlinear  internal  wave  motions.  Following  Bretherton 
and  Garrett  (19),  we  define  the  displacement 


£(r,t)  = Y(r,t)  - r 


[2.6] 
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and  consider  |c|  to  be  a small  quantity.  The  Lagrange  multiplier  P 
may  then  be  expanded  as 


P(Y,t)  = P(r,t)  + U*v)P(r,t)  + -J-  (5^7^  )P(r,t)  + ... 


[2.7] 


Substitution  into  [2.2]  yields 

L = / 7^-  [ 4-  P I C| 2 + Pg-C  - \ Pf*(C  * t) 

L0  ' - £ - * 

- U-V)P(r,t)  - (ti^j7i7j)P(r’t)  + •••  ] . [2.8] 

where  we  have  used  the  relation 


/ d3r  P(Y,t)J(3Y/3r)  = /d3Y  P(Y,t)  = / d3r  P(r,t)  , 


and  omitted  terms  independent  of  {. 

Next,  we  define  the  pressure  fluctuation 

ir(r,t)  = P(r,t)  - P(r)  [2.9] 

and  assume  it  is  a small  quantity  of  the  same  order  as  |c|.  On  intro- 
ducing [2.9]  into  [2.8]  and  collecting  terms,  we  obtain 
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d3r 


Li  = f r tp«*5  " (C-V)P(r) ] , 

X ^0  

L2  = / ^ [ jr  p\i\2  - \ f*U  x O - (c-v)t.  - \ (cicJvi7J)p(r)] 
l3  - / T7  1 - T < - T "iijWi'kWr"  • 


[2.10] 


The  equations  of  motion  in  each  order  are  obtained  by  variation  of 


rj = Lj dt 


with  respect  to  £ &nd  n.  Variation  of  1^  with  respect  to  £ gives 


pg 


= VP(r ) . 


[2.11] 


Thus  P is  a function  only  of  z and  is  identified  as  the  hydrostatic 
pressure. 

We  may  use  [2.1l]  in  the  Lagrangian  to  eliminate  P and  ignore  I 
from  now  on.  Then 


L = L2  + • • • » 


L = / ^ t T-  pUI2  - -5-  Pf-(C  x o - (C-V)ir  - p N 2d)  , 


L3  = /Tf  ^TWiV’-T^Ws1 


[2.12] 
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vhere  N2  is  defined  by  [2.1]. 

Linear  equations  of  motion  are  obtained  from  variation  of  I 
pi  + pf  x l + V*  + p n2C3z  = 0 , 

= 0 . [2.13] 

ft 

We  now  neglect  the  horizontal  components  of  f . This  enables  us 
to  separate  the  horizontal  and  vertical  parts  of  [2.13]  and  expand 
C3(x,z,t)  in  the  rectangular  area  zQ  as+ 

hu,z,t)  = f i V‘>  V 

a=l  k . 

A-k«  = J 

(the  first  sum  extends  over  all  positive  integers  o).  Our  task  will  be  to 
obtain  an  expression  for  the  Lagrangian  [2.12]  in  terms  of  the  field  ampli- 
tudes A^ft). 

Using  the  definition  [ 2 . lU ] and  some  straightforward  algebra, 
equations  [2.13]  yield 


•See,  for  example.  Ref.  (l),  p.  239- 

+Henceforth  al]  vectors  with  the  exception  of  f will  be  two-dimensional 
in  the  horizontal  plane. 
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V3pV3Wka  + pk 


N2-U2(k) 

u>2(k)-f2 

a 


W.  = ° , 

ka 


u>  (k)  > 0 , 

a 


\a<-K)  ' V°>  ' 0 • 


This  is  a Sturm-Liouville  equation  for  the  modefunctions  W,  (z)  and 

ka 

eigenvalues  wa(k).  The  orthogonality  relations 


£ lH  ' S«B 


are  readily  deduced  from  [2.15].  The  quantity  pQ  may,  for  example, 
chosen  as  p(0).  The  Fourier  amplitudes  satisfy  the  equation 


\«  * \a  ■ 0 


The  horizontal  displacement  is  expressed  in  the  form 


5h(*.z»t)  = I I [ik\a  + ~2  ^ka  6 


ik*x 


a=l  k 


where  §ka  = A^a  and  W'  = dW/dz.  The  pressure  fluctuation  is  then 


' - - » M <**„  * f2V  7?  "k« * 

a=l  k k 


ik*x 


[2.15] 


[2.16] 

be 


[2.17] 


[2.18] 


[2.19] 


Expansions  [2.lU],  [2.l8]  and  [2.19]  may  be  introduced  into  the 
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expression  [2.10]  for  to  give 


■s  ’ f I l 


a=l  k 


tf2(]t)_f2  (AkaA-ka  ' h’a(k)AkaA-ka) 


[2.20] 


To  obtain  the  nonlinear  equations,  we  first  write 


§ = e ^ + c <■  2 + . . . , 


IT  * E 7T^  + e *2  + • • • » 


[2.21] 


where  are  the  fields  given  in  the  linear  approximation  by  equations 

[2. lU ] through  [2.19],  and  e is  a small  parameter.  Now,  because  I 2 is 
stationary  for  small  variations  about  the  solutions  * , we  have 

12  = / L2  = e2  / L 2 (^1,ir1)dt  + 

13  = / 1*3  U.iOdt  = e3  / L3  ( ,1Ti )dt  + 0(eS.  [2.22] 

Thus,  since  we  do  not  wish  to  obtain  equations  valid  beyond  the  order  of 

3 

c we  may  use  the  linear  expressions  [2.lU]-[2.19]  in  deriving  L 2 and  L„. 

In  particular,  we  may  use  [2.17]  to  eliminate  time  derivative  terms 
from  L^.  It  is  important  to  realize  that  the  nonlinear  fields  will  con- 
tinue to  be  expressed  in  the  form  [2.lU ]— [2.19]  except  that  [2.17]  will 
be  modified. 

We  wish  to  express  our  equations  in  Hamiltonian  form.  As  described 
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in  some  detail  by  Meiss  and  Watson  this  can  be  conveniently  done  by 
replacing  the  Fourier  amplitudes  by  canonical  action  angle  variables 

^ka  * 0ka ' 


^ka  r 


ka 


\ 2u  e. 
a ka 


VaAT 

a ka 


0 


n — i0ka  rz i0-ka 

VJ ka  e ' - \J_V  e 


-ka 


-ie 


ka 


6 - " + V^ka  6 


ie 


■ka 


'ka 


2 ' 
0)  -f 
a 


[2.23] 


The  resulting  Hamiltonian  is  obtained  after  straightforward  but  tedious 
effort : 


H = H2  + H3  , 


H2  ■=  l «,<»>  Jka  . 

a ,k 


Vj  exp[i(ek-9t-0  )] 

k,£,m  ... 


+ 6V+*+m  r2(k’^ra)  exP[i(VV9m)1} 


+ c.c. 


[2.25] 


•See  ref.  (5),  PP-  296-323 


r 7 ^ 
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(The  coefficients  and  are  given  in  the  Appendix.) 

To  simplify  notation  in  we  have  written  k for  the  index  pair  (k,a), 
etc. 

The  equations  governing  the  evolution  of  the  internal  wave  field  are 
Hamilton's  equations: 


j - _ 9H 
ka  " “ 39 


ka 


: _ 3HL 

3ka  " 3J 


[2.26] 


ka 


For  linear  waves,  corresponding  to  neglect  of  H^,  we  have 


Jka  = ° » Ka  = Ua(k)  * 


[2.27] 


The  term  in  H describes  the  effect  of  all  interacting  triadB 
of  waves  which  satisfy  the  wavenumber  conservation  restrictions  k-t-m  = 0 
or  k+fc+m  = 0,  and  allows  for  transfer  of  energy  between  members  of  the 
triad.  Since  is  assumed  small  compared  with  Hg,  equation  [2.27] 
indicates  that  the  terms  involving  ^ are  rapidly  oscillating  on  the 
time  scale  of  energy  transfer  among  the  modes.  We  shall  henceforth 
neglect  these  terms,  setting  - 0 in  the  numerical  experiments  to  be 
described  later  in  this  paper.  A similar  argument  applied  to  the 
terms  suggests  that  only  those  wave  triads  for  which 


A = id  (k)  - u>„(0  - to  (m)  2r  0 
any 


[2.28] 


will  have  dynamical  significance  in  our  studies. 


r 


-13- 


We  shall  refer  to  [2.28]  as  a "resonance  condition."  Those  terms 
in  for  which  the  resonance  condition  is  not  met  can  be  transformed 
away  by  a sequence  of  canonical  transformations.  The  resulting  Hamiltonian 
has  the  form 


H 


1 


where  is  a function  of  the  transformed  action  variables  and  contains 
only  wave-wave  interactions  for  which  the  resonance  condition  [2.28]  obtains. 

For  our  present  purposes  we  shall  assume  that  those  triads  not  satis- 
fying the  resonance  condition  may  be  neglected.  How  large  the  resonance 
mismatch  A must  be  for  a triad  to  be  neglected  is,  of  course,  a question 
that  must  be  answered  quantitatively.  We  shall  return  to  this  point  in 
Section  I4. 

Another  set  of  convenient  variables  are  the  dimensionless  action- 
amplitude  variables^ 


[2.29] 


They  are  related  to  the  original  Fourier  coefficients  by 


. - ■ V!  , 

\a  1 2k  ~ 


-ka 


[2.30] 


+Nq  and  B are  constant  quantities  having  dimensions  of  frequency  and 
length  respectively.  They  may  be  considered  as  scale  parameters  of 
the  Vaisala  profile  N(z). 


r 
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In  the  linear  approximation,  represents  the  amplitude  of  a 
traveling  plane  wave  having  wave  number  k and  mode  number  a.  It  is 
normalized  to  describe  the  wave  slope  amplitude. 

The  equations  of  motion  in  terms  of  the  a.  are 


V = (aka*H) 


= - 2i 


3H 


3akc 


[2.31] 


Here  (a,H)  represents  the  Poisson  bracket  of  a with  H.  The  explicit 
form  of  equation  [2.31]  as  obtained  from  [2.25]  is 


\ + iwa\  = J 


dm 


dm 


6.  . G. 

k-i-m  dm 


a.  a ] 
d m 


[2.32] 


The  coefficients  G are  expressed  in  the  Appendix  as  functions  of  r^. 

In  this  equation  there  are  two  types  of  interaction  terms.  These  are 
denoted  sum  and  difference  interactions  according  to  the  wavenumber  con- 
servation relations  k=m±d  respectively. 
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3.  Test  Wave  Model 


To  provide  a simple  model  for  preliminary  study  we  introduce 
the  test  wave  system.  Here  a single  (k,a)  mode,  identified  as  the  test 
wave,  is  permitted  to  interact  with  an  ambient  field  consisting  of  the 
remaining  modes.  However  the  ambient  modes  are  not  allowed  to  interact 
among  themselves.  Test  wave  models  are  often  employed  for  calculating 
relaxation  rates  in  transport  theory,  where  standard  approximation 
schemes  in  statistical  mechanics  may  be  applied.  Numerical  integration 
of  equations  such  as  [2.26]  can  be  used  to  evaluate  the  validity  of  these 
approximation  methods.  Later  papers  in  this  series  will  study  this  problem 
in  detail. 

We  can  also  use  the  test  wave  model  to  test  whether  a particular 
background  wave  amplitude  spectrum  is  an  equilibrium  spectrum  with  respect 
to  the  dynamics.  If  it  is,  the  mean  action  of  the  test  wave  (averaged 
over  ambient  initial  conditions)  must  decay  to  its  equilibrium  value. 

Our  system  consists  of  2M+1  waves  formed  in  M triads.  We  let  k 
designate  the  test  wave  mode  (k,a)  and  !,  ,m  designate  ambient  wave  pairs 
(£,8),  (m,y)  respectively.  The  appropriate  equations  of  motion  are  then 


TEST  WAVE:  a,  + iu  a.  * J [«.  , G*  a a + * Gkta»a  ] 

K a K L k-£-m  tn  t n k+£-m  m i.  m 

i,m  _ _ „ _ _ _ 


BACKGROUND:  4 + ia>  a = 6„  . 2G*  a.  a + 6„  w G^\  a* 

l g l £-k-m  km  k m £-k+m  k K - 


+ 6 . 
£+k-m  m k m 


m 


[3.2] 
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It  is  easily  verified  that  [3.2]  admits  M+l  conservation  lavs.* 

In  our  computations,  initial  conditions  for  the  ambient  modes 
are  chosen  from  a Gaussian  distribution  for  the  action-amplitudes. 


?(\a>  ' 2^~exp  [-|*ka|2/2oka]  ‘ 

- *a 

Thus,  the  initial  wave  phases  are  f»~ > 1 1 *rtriPution  uniform 

in  -w,».  Averages  over  ambient  initial  conditions  with  distribution 

[3.3]  will  be  denoted  <>  . Itoe  mean  square  action  amplitude  o is 

ka 

related  to  the  pover  density  spectrum  for  the  vertical  displacement 

> = [ /d2k  t|/(k,a ) ; 

a=l 


< 


> 


Equivalently,  the  mean  wave  action  (see  equation  [2.29])  is  given  by 


<J 


ka 


> 


<*>f(k)-f2 


^(k,a)  . 


The  specific  iji(k,a)  used  for  this  study  was  the  Garrett-Munk  ”75+" 
spectrum  (17),  an  experimentally  derived  spectrum  for  internal  waves 
which  appears  to  be  an  equilibrium  spectrum  for  the  ocean. 


[3.3] 


[3-M 


[3.5] 


[3.6] 


•See  J.  D.  Meiss  and  K.  M.  Watson,  ref.  (5),  pp.  296-323. 
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^ . Numerical  Methods  and  Results 

Our  emphasis  in  this  section  will  be  on  computational  techniques 
and  qualitative  features  of  solutions.  Detailed  application  to  the  ocean 
will  be  published  separately. 

We  represent  the  continuous  wave  system  with  a finite  number  (2M+1) 
of  waves.  Each  wave  must  correspond  to  a wavelength  which  "fits”  into 
the  (assumed)  square  area  This  implies  a base  grid  in  wavenumber 

space  with  grid  spacing 


Ak  = Ak 
x y 


2n 


[U.l] 


Since  the  test  wave  wavenumber  k is  fixed  (given),  the  wavenumber 
conservation  equations 

k = ml*,  [U.2] 

imply  that  a single  grid  can  act  as  coordinate  system.  For  example, 

chossing  the  x-axis  along  k,  each  point  (m  ,m  ) defines  a triad  as  shown 

- x y 

in  Figure  1.  As  discussed  in  Section  2,  the  most  important  triads  for 

energy  transfer  are  expected  to  be  those  for  which  the  resonance  condition 

A=oj  + u>„  - u>  -0  holds.  For  a given  set  of  modenumbers  a,6>Y 

a 8 Y 

and  a given  resonance  mismatch  A,  each  resonance  condition  defines  a 

curve  in  the  m ,m  plane.  Examples  of  curves  for  various  A are  shown 
x y 

in  Figure  1.  Our  procedure  will  be  to  include  those  triads  which  lie 

on  grid  points  for  which  |a|sA  • For  a given  computation  the  parameters 

max 
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Ak  and  Amnv  must  be  obtained  empirically.  Before  discussing  this  further 
we  present  a sample  computation  of  the  full  system  of  equation  to  illustrate 
the  general  character  of  the  solution. 

! 

Figure  2 displays  a typical  example  of  the  mean  test  wave  action  for 

the  ocean  internal  wave  system.  The  frequency  and  modenumber  of  the  test 

wave  are  u =7f  and  a=U,  respectively.  We  include  the  first  nine  WKB 
a 

vertical  modes  of  the  ambient  spectrum.  This  corresponds  to  including 
fifteen  resonance  curves  in  the  calculation.  The  magnitude  of  the  initial 
test  wave  action-amplitude  is  chosen  to  be  0.1  ^ and  the  corresponding 

initial  phase  is  fixed  arbitrarily.  The  ambient  initial  conditions  are 
picked  with  a random  number  generator  according  to  the  distribution  [3.*»]. 

To  integrate  the  equation  of  motion  [3.2],  we  use  a variable  step  size, 
variable  order  Adams  integration  scheme  (20).  The  test  wave  action  is 
averaged  over  1=100  ambient  realizations,  integrating  to  t=1000/NQ  for 
each.  The  entire  computation  required  approximately  7 minutes  C.P.  time 
on  a CDC  7600  computer.  The  salient  features  of  Figure  2 are  a rapid  rise 
of  the  test  wave  action  by  a factor  of  100  in  the  first  20  time  units, 
followed  by  small  oscillations  about  a value  slightly  larger  them  its 
expected  (Garrett-Munk)  equilibrium  value  (indicated  by  the  dotted  hori- 
zontal  line).  We  have  found  that  these  fluctuations  decrease  with  increasing 
I,  emd  that  the  mean  value  is  constant  upon  variation  of  the  initial  test 
wave  amplitude  providing  this  amplitude  is  small  (<0.25  ^ * 

The  effect  of  choosing  different  values  of  Ak  and  Amax  is  best  illustrated 
by  a simple  model  calculation.  For  the  remainder  of  this  section  we  keep  only 
the  sum  interactions  in  [3.2]  and  arbitrarily  set  the  coupling  coefficients 


w 
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to  = G^m  = G™*'  » G = -4-  (-1+i),  ignoring  the  complicated  wavenumber 

x.m  K k d 

and  frequency  dependence  of  the  read  ocean. 

To  investigate  the  effect  of  off-resonant  triads  ve  consider  25 
triads,  13  with  A=0  and  12  with  A=A  _•  Figures  3 and  *4  present  plots 

IQEX 

of  the  short  time  behavior  of  the  mean  action  for  various  A . In  Figure  3 

max 

the  phases  of  the  ambient  modes  at  t=0  are  chosen  so  that  the  initial 

growth  rate  of  the  test  wave  action  is  maximal.  That  is,  the  sum  of  the 

phases  of  the  background  waves  in  [3.2]  plus  the  phase  of  G is  zero.  It 

is  known  that  this  set  of  initial  conditions  results  in  "square  wave" 

phase  oscillations  and  maximum  coupling  of  the  interacting  waves  (21). 

These  coherent  solutions  in  Figure  3 show  that  the  effect  of  off-resonant 

triads  is  to  produce  oscillations  about  the  growth  rates  of  the  resonant 

triads.  For  t<<l/AmHy  the  growth  rate  is  essentially  given  by  the 

A =0  curve . 
max 

For  later  times  the  action  oscillates  about  the  A =“>  curve  (where 

max 

in  effect  only  13  resonant  triads  are  kept,  all  with  A=0).  The  frequency 
of  these  oscillations  depends  linearly  on  A 

max 

In  Figure  U the  initial  phases  are  chosen  randomly.  From  [3.2]  we 

see  that  the  intial  growth  rate  of  the  mean  action  is  zero.  For  later 

times,  the  growth  rate  (calculated  as  an  average  over  100  ambient  initial 

conditions)  is  about  one  half  of  the  maximal  rate.  As  A is  increased, 

max 

the  growth  rate  decreases  as  in  Figure  3,  but  now  the  oscillations  are 
washed  out  by  the  phase  averages.  The  ambient  phases  do  not  remain  com- 
pletely random  for  t>0:  If  we  make  the  random  phase  approximation  using 
the  Krylov-Bogoliubov-Mitropolsky  (1*4 ) (KBM)  perturbation  theory  to  obtain 
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an  estimate  for  the  growth  rate,  we  get  a rate  that  is  less  than  that  of  the 
dynamical  calculation.  The  KBM  and  maximal  rates  represent  extremes, 
the  first  implying  complete  incoherence  while  the  latter  implies  complete 
coherence.  We  plan  to  quantify  the  significance  of  phase  coherence  in  a 
later  paper. 

As  a final  model,  more  closely  related  to  the  full  ocean  calcula- 
tion, we  include  only  triads  with  modenumbers  a=1»,  6*6  and  y=2  corres- 
ponding to  the  resonance  curves  of  Figure  1.  Again  we  choose  the  G's 
to  be  equal.  With  this  model  we  can  discuss  the  effect  of  both  Ak  and 
Amax  °n  the  rate.  For  this  purpose  we  tentatively  assume  that  the 

test  wave  relaxes  exponentially  to  equilibrium.  This  type  of  relaxation 
would  be  predicted  by  a Langevin  model  for  the  system.  To  provide  a rough 
estimate  for  the  growth  rate,  v,  we  calculate  the  time,  t^  , required  for 
the  test  wave  to  grow  to  half  of  its  equilibrium  value.  Then 
v = -log2/t^  . The  exponential  relaxation  assumption  is  clearly  not 

2 2 

valid  for  very  short  times  when  <'|ajeal  ^ grows  as  t (since  the  ambient 
phases  are  initially  random) . For  longer  times  the  oscillations  caused 
by  non-resonant  triads  must  also  provide  an  error  in  the  fit.  However 
for  the  purpose  of  investigating  the  effects  of  changing  Ak  and  A 

max 

our  measure  is  sufficient. 

In  Figure  5 we  plot  v as  a function  of  A for  various  values  of 

max 

Ak.  These  curves  approach  v » .02Nq  for  Amax  £ »1Nq  independently  of 
Ak.  Since  varying  Ak  changes  the  number  of  triads  within  a given  reso- 
nance width,  we  see  that  keeping  a small  number  of  triads  is  sufficient 

to  obtain  a good  estimate  of  v.  For  example,  the  value  of  v for  A *.l 

max 
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and  Ak=.27  (ll*  triads)  differs  by  only  15%  from  that  with  Ak=.lO 

. _1* 

(93  triads).  Varying  the  initial  action  between  10  o and  .050. 

ka  kot 

does  not  significantly  change  the  growth  rate  or  the  final  action. 
Further  increase  of  the  initial  action  changes  the  growth  rates  drasti- 
cally indicating  that  the  exponential  assumption  breaks  down. 

Figure  5 indicates  that  quantitative  results  may  be  obtained  for 

A sr  5v.  Since  A and  v define  the  only  available  time  scales, 
max  max 

these  quantities  can  be  expected  to  scale  together  and  produce  reliable 

results  if  A is  chosen,  self  consistently  to  be  "a  few  times  v." 
max 

Finally,  in  Figure  6 we  see  that  the  time  required  for  the  test 
wave  action  to  reach  equilibrium  does  not  depend  strongly  on  A 

max 

Although  initially  the  two  curves  grow  at  quite  different  rates 

[v(A  =.03Nr.)  *=  1.5  v(A  =.2Nn)],  for  later  times  the  curves  grow 
max  u max  u 

at  similar  rates  and  reach  their  equilibrium  values  at  comparable  times. 
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The  coupling  coefficients  which  appear  in  equations  [2.25]  and 
[2.32]  of  the  main  text  are  defined  explicitly  by  the  following  set 
of  equations: 


f2  = Pq(  32u)a<»)giD^ekae|lgem^ ) ^ [G(k;i,  ,m)  + G(t.;m,k)  + G(m;k,i.  / ]/3, 


[A.l  ] 


^ « xv  ^ 

G(k;!,m)  = (u>2-w2)  Ixm  • — ( — ~ - 0^2*0 

B y - - “a  \ k2  k2tm  ' ' 


(«•* 

,V.2  \ - - 


k)(k*m)  + (fc*m)‘ 


[A. 2] 


r,  (k,(u  ; l ,<dq  ,m,w  ) = 3r9(k,u)  ) , 

1 . a . p . Y d ~ 0 “ T 


[A. 3] 


^ - -i—  fdzoW'WW  v = f dz  p W'  W'  W'  , 
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Figure  1.  Resonance  curves  for  sum  interaction  using  the  deep 

2 2 2 
ocean  dispersion  relation:  w^vk)  = f + [kBNg/ii(a-Ji)  ] . 

For  these  curves  a = U,  B = 6,  y = 2,  ^(k)  = 7f. 

Arrowed  is  a typical  wavenumber  triad. 


Mean  test  wave  action  as  a function  of  time  with  91  triads 
on  15  resonance  curves  and  coupling  coefficients  for  the 
deep  ocean,  a = k,  u>a  = 7f  and  the  mean  is  taken  over 
I = 100  initial  conditions. 

Mean  test  wave  action  for  the  two-row  model.  The  ambient 
initial  phases  are  fixed  to  give  maximum  initial  growth 
rates.  I = 100. 


Figure  U. 


Figure  5- 


Figure  6. 


Mean  test  wave  action  for  the  two-row  model.  Ambient 
initial  phases  are  chosen  randomly  for  the  dashed  curves. 
Curves  labeled  "maximal"  are  from  Figure  3.  The  KBM 
curve  assumes  the  random  phase  approximation. 

Growth  rate,  v,  as  a function  of  A and  Ak  for  the 
, infix 

single  resonance  curve  model. 

Mean  test  wave  action  for  the  single  resonance  model.  For 

t < 100/N-  the  A solution  grows  faster  than  the  A = .03 
u in  8.x  max 

solution.  For  longer  times  (right  scale)  evolution  is  at  a 


similar  rate. 
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